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ABSTRACT 



Context. The wings of the Can H and K lines provide excellent photospheric temperature diagnostics. At the Swedish 1 -meter Solar 
Telescope (SST), the blue wing of Can H is scanned with a narrowband interference filter mounted on a rotation stage. This provides 
up to O.'IO spatial resolution filtergrams at high cadence that are concurrent with other diagnostics at longer wavelengths. 
Aims. The aim is to develop observational techniques that provide photospheric temperature stratification at the highest spatial reso- 
lution possible and use them to compare simulations and observations at diff'erent heights. 

Methods. We use filtergrams in the Can H blue wing that were obtained with a tillable interference filter at the SST. Synthetic obser- 
vations are produced from three-dimensional (3D) hydro and magneto-hydrodynamic numerical simulations and degraded to match 
the observations. The temperature structure obtained from applying the method to the synthetic data is compared with the known 
structure in the simulated atmospheres and with observations of an active region. Cross-correlation techniques using restored non- 
simultaneous continuum images are used to reduce high- altitude, small-scale seeing signal introduced from the non- simultaneity of 
the frames when diff'erentiating data. 

Results. Temperature extraction using high-resolution filtergrams in the Ca ii H blue wing works reasonably well when tested with 
simulated 3D atmospheres. The cross-correlation technique successfully compensates for the problem of small-scale seeing diff'er- 
ences and provides a measure of the spurious signal from this source in diff'erentiated data. Synthesized data from the simulated 
atmospheres (including pores) match well the observations morphologically at diff'erent observed heights and in vertical temperature 
gradients. 



Key words. Techniques: image processing - Sun: faculae, plages - Sun: granulation - 
Techniques: high angular resolution 



Sun: photosphere - Sun: magnetic topology - 



1. Introduction 

While the Ca ii H and K lines have been used for chromospheric 
studies for a long time, their usefulness as a photospheric tem- 
perature diagnostic remains largely unexplored in filtergrams 
with present-day large- aperture telescopes. 

The extended damping wings of Ca ii H and K are insensitive 
to line-of- sight velocities and magnetic fields. They provide in- 
formation over a wide depth range, are formed in Local Thermal 
Equilibrium (LTE), and allow data to be recorded at close to 
the highest spatial resolution achievable with ground-based op- 
tical telescopes. The absence of Doppler signals is an advantage 
when extracting temperature information. Furthermore, the up- 
per photosphere processes that are sampled in the inner wings 
of these lines provide important clues to studying the chromo- 
sphere. Beck et al. ( 2008) made a good case for this and used 
spectra in the inner Ca ii H blue wing to make a comparison with 
theoretical heating models and to identify photospheric waves as 
a driver for chromospheric line-core emission. 

Slit spectroscopy in the wings of these lines has been used 
to successfully constrain semi-empi rical one-dimensional model 
atmos pheres f or the "average Sun" dAyresll975l:lAyres & LinskyI 
11976'; Waltori'1987') as well as for small magnetic feature mod- 
els ( Fang et al. 2001)^ Also targeting bright magnetic features, 
IShine & LinskyI (Il974l) used the Can K line wings, the fact 
that Can is mostly in the ground state throughout the forma- 



tion height of the line, and the assumption of hydrostatic equi- 
librium to derive the temperature as a function of column mass 
from observations. Th is method was used for late-type stars by 
lAvres & Linskvl (1 19751) and Avr es et al] (1 19741). It was als o used 
for a sunspot penumbra by Rouppe van der Voori (l2002h . who 
investigated its temperature structure using a scan obtained with 
the slit spectrograph at the former 50-cm Swedish Vacuum Solar 
Telescope, thus obtaining a 3D temperature model of the penum- 
bra approaching ^/23 resolution. That sa me data and method 
were again used by lSheminova et al.] (12005) to compare the com- 
puted stratification obtained at diflferent bright points, with pre- 
dicti ons made for vert ical two-dimensional thin fluxtube mod- 
els. Sheminoval (1201 2h produced synthetic spectra from hydro- 
dynamic (HD) and magneto-hydrodynamic (MHD) simulations 
and investigated the accuracy of a "manual" inversion technique 
for extracting temperature stratification, demonstrating good re- 
sults for selected small structures. 

Recent high-resolution work based on filtergrams recorded 
in the wings of these lines has focused on qualita- 
tive interpretation of images in th e Can H blue wing. 
iLeenaarts & Wedemeyer-Bohml (l2005l) observed the inner blue 
wing of Can H (396.667 nm) an d co mpared it with l i ne cor e fil- 
tergrams. Scharme r et a D (I2007h and'Langhans et all (l2007l) also 
used an inner wing position to ex plore the upper photosphere of a 
sunspot penumbra. iBerger et al.l ([2004b observed in the G-band 
and close to the Ca ii H line core to characterize bright points 
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and small pores at (7/10 resolution. This last work used a 0.29 
nm wide filter in a setup at the Swedish 1 -meter Solar Telescope 
(SST) that is similar to the one used in the present paper regard- 
ing the secondary cameras. 

Few high-resolution semi- empirical 3D models o f solar fine 
structure have been published. ISocas-Navarrol (1201 ll) . recogniz- 
ing the lack and the need for these in the context of abun- 
dance determination, used Hinode spectra in the 6301 and 6302 
line region with th e non- LTE (NLTE) inversion method of 
ISocas-Navarro et aP d 19981) to produce a model of quiet gran- 
ulation at disk center at 0'.'3 resolution. 

At the SST, O'.'IO is achievable at the wavelength of the 
Can H and K lines, and the use of filtergrams allows for control 
and compensation of seeing eff'ects using imag e reconstruction 
techniques. In this work, we apply the method o fiShine & LinskyI 
(11974 to filter grams in the blue wins of Can H. We address 
the main image processing issues with such an adaptation of 
their method to filtergram data and test the method by producing 
and analyzing synthetic observations from an HD simulation. 
Finally, we compare observations and simulations and discuss 
diff'erent photospheric structures in diff'erent heights and in ver- 
tical temperature gradients. 

2. Setup and observations 

The observations were perf ormed on 23 May 2010 with the 
SST (IScharmer et al.ll2003a|). including the 3 7 - electrode-mirror 
adaptive-optics system (" Scharmer et al.ll2003bl) . The target was 
active region NOAA 11072, including a sunspot, pores, and as- 
sociated bright points (see Fig.[T]). The field of view is 62'' x 62'' 
with an image scale of 0'.'034 per pixel (25.0 km per pixel on the 
Sun). The active region was observed at a heliocentric angle of 
~15°corresponding to ji = cosO = 0.97. 

The setup, known as the blue tower, included four Megaplus 
II es4020 cameras. The camera fed by the first beamsplitter 
cube, thus receiving the most light of the four cameras, was 
placed behind a 0.11 nm interference filter mounted on a rota- 
tion stage. The central wavelength at normal incidence of this 
filter is 396.88 nm, just red- wise of the core of Can H, but ro- 
tating the filter shifts the center of the pass band towards the 
blue. This provides a simple way to scan the entire blue wing 
of the Ca ii H line at the cost of lower spectral resolution and 
image degradation for large tilt angles, as explained by Beckers 
(1998). The second camera was mounted behind a 0.12 nm wide 
filter at 396.47 nm in the Can H blue wing. The wing camera is 
included in the M ulti-Object Multi-Frame Blind Deconvolution 
(Ivan Noort et al. l2005. MOMFBD) reconstructions but not in the 
present analysis because we focus on the tillable filter. The last 
two cameras were mounted behind a wideband (WB) filter with a 
FWHM of 1 nm centered on 395.37 nm in the pseudo continuum 
between the Can H and K lines. One of these was intentionally 
defocused by 7 mm for approximately one wave peak to peak of 
phase diversity (PD). 

The present data are from a single scan, which consists of 
91 images per camera with the tunable filter camera scanning 
seven wavelengths in the Can H blue wing plus line core (13 
images per tunable filter wavelength). The scanned wavelengths 
and diff'erent filter-position properties are listed in Table [T] The 
exposure time was 8 ms, with the full scan taking 9.7 seconds. 

The SST setup allows for co-observati ons in the red with 
the CRISP Imaging Spectropolarimeter ( S charmer et al.l [20081 
CRISP). The data, acquired at 14:11 UT, are nearly simultane- 
ous with C 1 538.0 nm, Fe 1 630.2 nm, and 630. 1 n m observations. 
Analyses of these observations are published in IScharmer et al.l 



(1201 ll) and IScharmer & HenriquesI (|2012|) . where convective 
downflows were detected in the sunspots penumbra. We intend 
to combine the Ca ii H and CRISP data sets in future work. 

3. Data reduction 

Data reduction is centered around i mage reconstruction using 
the MOMFBD code ( van Noort et al.>2005) with PD. MOMFBD 
was used with the 36 most significant Karhunen-Loeve coeffi- 
cients. Since we aim at diff'erentiating data, additional steps were 
taken to properly account for the complications that arise from 
the tillable filter and to address frame non- simultaneity. The data 
reduction procedures are summarized below in chronological or- 
der and then discussed in detail in the specified sub-sections: 

- Non-linearity correction for one of the cameras (Sect. 13.11) 

- Camera alignment with deconvolved pinholes (Sect. 13.21) 

- Dark correction and flat fielding 

- MOMFBD with extra image objects (Sect. 13.31) 

- Compensation for tilted-filter PSF eff'ects (Sect. 13.21) 

- Dewarping using ext ra im age objects (Sect. 13.31) 

- Noise filtering (Sect. 13.21) 

3.1. Non-linearity correction and ghosts 

The MegaPlus II es4020 camera CCDs consist of two halves 
with separate read-out circuits. One of the halves of the PD 
camera showed a consistent non-linearity in the gain. After 
MOMFBD reconstruction, this diff'erence between the two 
halves gave rise to a strong central artifact. This non-linear re- 
sponse was measured by Dr. P. Siitterlin (2010 private communi- 
cation) and was well modeled by a six-degree polynomial. The 
non-linear gain would introduce, at its peak, ~ 6% extra counts 
compared to the linear response curve from the other half. For 
this data set, that polynomial was used to correct each individual 
raw science and flat-field frame after dark correction for the af- 
fected camera. The camera was repaired between the 2010 and 
201 1 campaigns and no longer presents this issue. 

In pinhole array images, we detected defocused ghost images 
of around 1% for all positions of the tunable filter. Based on the 
progression of their distance to the main image with increasing 
tilt angle, these likely come from a secondary image generated 
by inner reflections in the tilted filter. We ignore these ghosts. 
They will have an impact similar to a small amount of scattered 
light. Extreme inverse filtering and inspection of the images in 
high contrast in the umbra do not reveal any visible ghosts in the 
final reconstructed images. 

3.2. Compensating for tilted-filter effects 

The degradation eff'ects from til ting an inter ference filter in this 
setup are discussed in Lofdahl et al.l (1201 lb . where a compen- 
sation scheme in the context of MOMFBD reconstruction was 
also established. Here we apply that post-MOMFBD compensa- 
tion scheme, which includes deconvolution of the reconstructed 
images with the known elongated point spread functions (PSFs) 
from tilting and extra noise filtering. This compensation is pos- 
sible since the eff'ects of a tilted filter are separable from phase 
aberrations and thus do not inte rfere with wavefront sensing in 
MOMFBD (lLofdahletal.1120111) . We follow the pre-MOMFBD 
pinhole deconvolution as proposed in that work and addition- 
ally apply a noise filter to the deconvolved pinhole images in 
the same way as for each tunable-filter science position prior to 
MOMFBD. 
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Fig. 1. Wideband image of the whole field of view. Each tick mark is 5'' 
corners are used for intensity calibration. 



The highlighted "quiet Sun" areas in the upper left and lower right 



For the individual temperature maps, the images at each 
tilt-angle position were low-pass filtered with a diflTerent, 
individually optimized spatial filter, which was estimated from 
the noise levels of the final deconvolved images. However, when 
diflTerentiating temperature maps (as done in Sect. O, the data 
to be included in the computation need to be filtered in the same 
way in order not to introduce artificial signal due to diflferent 
resolutions of the two temperature maps. Since higher tilt angles 
require a more aggressive noise filtering (a lower pass threshold 
in the Fourier domain), we separate the data in two sets before 
diflTerentiating data. For each set, the same noise filter is applied 
to all narrowband (NB) images. This filter is built by using 
the minimum, at each coordinate point in Fourier space, of the 
individually optimized noise filters. It is close to the noise filter 
from the highest tilt angle image used in each set. For one of the 
sets we use all the filter positions requiring aggressive filtering 
and limiting the achievable resolution to about (X'lT. For the 
other set, we use only the smallest four tilt angles so that the 
analysis of the gradients in the upper photosphere can be done 
at the maximum resolution (see Table [2] for the different sets). 



3.3. Making non-simultaneous frames consistent 

It is not possible to perfectly compensate for seeing using im- 
age reconstruction. The diflferent, non- simultaneous, tunable fil- 



Table 1. Filter properties at different positions. 



Pos. 


Tilt 




FWHM 




Mi 


IgTsoO 


(0 


angle 


(nm) 


(nm) 


(nm) 


(g cm-2) 




(1) 


1?2 


396.84 


0.12 


396.81 


0.17 


-2.62 


(2) 


2?5 


396.74 


0.14 


396.71 


0.90 


-1.31 


(3) 


3?0 


396.67 


0.16 


396.64 


1.40 


-0.96 


(4) 


3°1 


396.57 


0.18 


396.51 


2.18 


-0.74 


(5) 


4?2 


396.47 


0.19 


396.43 


2.61 


-0.42 


(6) 


4?8 


396.34 


0.22 


396.33 


3.04 


-0.37 


(7) 


6?4 


395.93 


0.28 


395.88 


3.56 


-0.27 



Notes. The first column lists the position number as used in the text. 
Xc is the central wavelength of the passband for each tilt angle and X^^i 
the effective wavelength defined in Sect. 14.11 mi is the average depth 
for the whole field of view computed using Eq. ID and Ig T500 is the 
corresponding optical depth value in HolMul. 



ter reconstructions will be affected by residual differential see- 
ing effects (in particular, small-scale geometric distortions) that 
translate to false signal when differentiating data. However, it is 
possible to improve consistency between the non- simultaneous 
NB objects by using a simultaneous constant reference. Such a 
reference is provided by the wideband object (shown in Fig.[3]in 
light green). 

In order to take advantage of this reference, the MOMFBD 
reconstruction was done with one additional object ("WB" in 
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Fig. 2. Different transmission profiles (Fi (A)) of the interference filter for the tilt angles used. The FTS atlas (lNeckeliri999h is shown in red, 
together with the synthetic spectrum from HolMul with blends (black continuous line) computed at ytz = 0.953. The filter profiles are normalized 
to the peak transmission of the narrowest profile in the plot, over the line core, that corresponds to the smallest tilt angle used. The large blend at 
396.15 nm is from All. 



Fig. O for each NB object consisting only of wideband frames 
that are simultaneous with individual frames for a particular NB 
object ("NB" in Fig. [3]). The resulting reconstruction thus pro- 
duces one WB image per NB reconstructed image that is affected 
by the same seeing and same reconstruction PSFs as the NB im- 
age. These individual WB images can then be used together with 
the reference (anchor) WB object spanning the full- scan time 
length to compute a geometrical distortion ("dewarp") matrix 
obtained by using cross-correlation techniques to measure small- 
scale diflTerential image motion. Applying that dewarp matrix to 
the corresponding NB object improves the consistency between 
that NB object and all the other non- simultaneous NB objects, 
even though the procedure does not improve small-scale differ- 
ential blurring. 

For this work the dewarping procedure involved applying 
three dewarp matrices in succession, with the last using 16 by 
16 pixel boxes ((y.'544) for correlation. This is approximately 5.5 
times the diffraction limit (Off 10) that we believe is achieved in 
the WB objects used as well as for filter positions (1), (2), (3), 
and (4) for this data set. Outliers in each dewarp matrix were re- 
moved, the matrix was smoothed, and a box overlap of one-third 
was used. Due to the high contrast of the images and the abun- 
dance of visible small structure, perhaps even smaller boxes for 
correlation could be used. The extra WB objects were given zero 
weight in the MOMFBD wavefront determination. 

To evaluate both the efficacy of this method, and provide a 
pixel-by-pixel measure of seeing-induced signal when diflTeren- 
tiating data, we compute a difference map between each of the 
time- segmented WB objects and the anchor WB object. Using 
these difference maps, we can assess the amount of signal in- 
troduced by non- simultaneity and estimate how much false sig- 
nal the dewarping procedure reduces. The full-field average im- 
provement is small, as would be expected in good seeing condi- 
tions. However, there are marked improvements in areas that are 
related to high horizontal intensity gradients. This is the case for 
the filament heads shown in Fig.|4l Around the central filament, 
the full range of the difference map before the dewarping proce- 
dure is -4.5% to 3.6% of the WB dynamic range and improves 
to -3.6% to 2.5% after dewarping. 

When computing temperature gradient maps, we used data 
compensated in this way. When looking at maps from single im- 




Fig. 3. MOMFBD scheme with the extra WB objects, exemplified for 
two tunable filter positions. For each NB object (marked as "NB"), a 
collection of simultaneous WB images (marked as "WB Extra 1" and 
"WB Extra 2") are selected and input as an extra object. The dark green 
images on top are from the WB PD camera. The scheme used for this 
work includes seven NB objects and their corresponding seven WB ex- 
tra objects, each with 13 realizations. 

ages, we apply no dewarping (such as in the intensity maps of 
Fig. |7] and Fig.©. 



4. Calibration and inversion method 

4.1. Intensity calibration 

The counts of each pixel were converted to absolute inten- 
sity using a synthetic profile computed in NLTE using the 
line- synthesis code MULTI (Carlsson 1986, version 2.3). The 
Holweger-Miiller ( Holweger & Mueller 1974, HolMul) model 
atmosphere is used because this well-tested empirical model is 
constructed to reproduce spectral lines in LTE. Bl ends were in- 
cluded in LTE from the VALD databa se (Piskuno v et al.lll995l: 
iRvabchi kova et al.lll997l: iKupka et al.l l2Q0a .1999). We also in- 
clude the Can K line as a background opacity for the full 
wavelength range. The obtained profile matches the National 
Solar Observatory Fourier Transform Spectrometer (FTS) atlas 
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Fig. 4. Left: WB image. Center: difference between an unwarped WB extra object and the WB anchor object. Right: same as center but using the 
dewarped sub-object. Each tickmark is 1". The scale in the center and right panels is the same. 



(lNeckellll999h reasonably well (see Fig. [2]) and provides infor- 
mation about the opacity contribution from Ca ii H and the back- 
ground sources at different wavelengths. Profiles were computed 
for a range of ji values, including that of the observations. For the 
observational quiet Sun reference, we selected two regions at op- 
posite extremes of the field of view as indicated in Fig.[T] These 
regions were chosen for their relatively low activity using an im- 
age from an inner wing filter position. An intensity calibration 
factor per filter position was then computed from these regions 
and the calibration synthetic spectra. For clarity, the expression 
for the calibrated observed intensity is 



Nr 



(1) 



where / is the filter position, Ni is the number of counts in the 
observed pixel, ^/,qs is the mean of the selected quiet Sun pixels, 
Fi (A) is the filter transmission profile ( see Fig.O computed from 
the tilt angle as in lLofdahl et al.l (1201 ll) . and (A) is the synthetic 
spectrum. 



4.2. Temperature and depth estimates 

Each filter position only provides one temperature and one wave- 
length point per pixel. The asymmetric intensity distribution of 
the line over each filter profile will cause the observed data to 
have effective central wavelengths that are blue- ward of the cen- 
tral wavelength values of each passband. For this reason, we 
multiply the filter profiles with the synthetic spectrum and use 
the wavelength of maximum intensity (/Im,/) when determining 
temperature, depth, and background opacities. This is verified to 
be at a wavelength not dominated by a blend. 

We compute the brightness temperature (Th) using the cali- 
brated intensity (//) and Aj^j- 



Nn(AA) 



5000 y 



0.39 



Trad 

Fw 



Trad 

Fw 



(5000) 



1/2- 



(3) 



(4) 



0.3155 and rad iative damp- 
\ both from iTheodosioul 



with an oscillator strength of fiu = 
ing constant of yrad = 1.42 x 10" s 
(1989). The coefficient Fw = 1.33 x 10"^ is calculated by 
IkouDpe van der Voort ( 2002) from the broadening cross-section 
of Barklem & O'Mara (1998). The neutral gas mean molecular 
weight is taken to be yUmoi = 1-29. Abundances, including that of 
Ca (Aca abo ve) and for th e elements used in MULTI, are taken 
from Asplund et al.l (l2009h . The continuum optical depth correc- 
tion (1 - Tc) was computed with MULTI using the background 
optical depth ratio to total opacity at T(/lm,i) = in the HolMul 
model. The core wavelength is taken to be Aq = 396.846 nm 
and AA = Aq - Aj^j- We use the one-p oint quadrature weight 

= 1.6 as defined and estimated by 'Rouppe van der VoortI 
(2002). Finally, mn, g, e, rrie, k, c are the hydrogen's mass, so- 
lar surface gravity, electron charge, electron mass, Boltzmann's 
constant, and speed of light, respectively. These equations as- 
sume LTE, and hydrostatic equilibrium and do not account for 
magnetic forces. 

The temperature dependence of the depth scale leads to each 
pixel being assigned a slightly different column mass. In order 
to compute temperature maps for a uniform depth-per- filter po- 
sition, the temperature stratification per pixel, as given by the 
seven images and the equations above, is interpolated with a cu- 
bic spline to a common depth scale. 



(2) 4.3. Temperature gradient maps 



where is the inverse Planck function. Then, assuming Tkin = 
Tb, we assign a column mass den sity for each pixel using ex- 
pressions (4) and (5) fro m Rouppe van der Voort ( 2002), who 
clarified the expression of lShine&Linskvl (11974 . and simply 
adapted to Can H: 



Each gradient map in Fig. [7] and Fig. [8] was computed from two 
temperature maps: 



Gij(x,y) = 



Tm,{x,y)-Tjnj{x,y) 



(5) 
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Table 2. Shown gradient maps, depths 



Number 


Using T maps 


m 


Resolution 




from positions 


(g cm-2) 




Gi2 


(l+),(2+) 


0.54 


a.^o 


G23 


(2+),(3+) 


1.15 


a.ao 


G34 


(3+),(4+) 


1.78 


a.^o 


G45 


(4-),(5-) 


2.39 




G56 


(5-),(6-) 


2.82 




Ge? 


(6-),(7-) 


3.30 




G57 


(5-),(7-) 


3.09 





Notes. The first column is an identifier for the gradient map as used in 
the text and figures. The second column identifies the filter positions 
used. The sign identifies images filtered with the low-resolution 
noise filter and the "+" sign images filtered with the high-resolution 
noise filter. The third column (m) is the depth in column mass. The last 
column is the maximum resolution allowed by the noise filter used. 



However, since the temperature map for each depth was in- 
terpolated to a single mean depth using a bicubic spline includ- 
ing the four nearest neighbors, gradients do include some in- 
formation from other sampled wavelengths and smooth out the 
quantities by a small amount. In Table [2] we show the depths cor- 
responding to the displayed temperature gradient maps and the 
temperature map pairs used to compute them. The images used 
to compute the temperature maps are filtered using two differ- 
ent sets of noise filters as described in Sect. 13. 21 and identified in 
TableEl 
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400 300 200 100 50 



400 300 200 100 50 

km 

Fig. 5. Temperature stratification for two columns, extracted from the 
synthetic images, are shown. The solid lines are the respective known 
temperature stratifications from the 3D HD simulated atmosphere. The 
dark crosshair marks the brightest near-continuum pixel and corre- 
sponds to the upper set of points, upper line, and the upper geometrical 
height scale. The white crosshair marks the darkest pixel and corre- 
sponds to the bottom set of points, bottom line, and bottom geometrical 
height scale, with zero defined by the mean optical depth (Igrsoo) = 1. 



4.4. Synthetic data 

Synthetic data were produced from two diflTerent simulations. To 
test the method and assist interpretation of non-magnetic areas in 
the observations, we use a 3D HD simulation snapshot computed 
with the code of IStein & N ordlund (1998) and kindly provided 
by Dr. M. Carlsson. To assist interpretation of areas containing 
numerous bright points and pores, we use a 3D MHD simu- 
lation snapshot from Lagerfjard & Nordlund (in prep.), kindly 
provided by Dr. A. Nordlund. This simulation has a mean total 
magnetic field of approximately 560 G around Ig T500 = for the 
section shown in Fig.[8l 

In each case, we again use MULTI to compute a synthetic 
profile with one plane-parallel atmosphere per pixel ("1.5D ge- 
ometry"), using a six-level Can model atom and the same pa- 
rameters (such as abundances, oscillator strength, and damping 
constant) used for the HolMul reference spectrum and Eq. ([3]) 
and including the same line blends. We take the same yu as the 
observations. Synthetic observations are then produced by mul- 
tiplying the synthetic line profiles with the diflTerent filter pro- 
files. The total area of the simulations is 6x6 Mm for the HD 
simulation and 12x12 Mm for the MHD simulation. Both snap- 
shots extend up to about 500 km above the mean continuum- 
forming layer of (Ig T500) = 1 with the HD simulation extending 
to 2.5 Mm below and the MHD simulation 20 Mm below the 
surface. Neither simulation has a chromosphere. The number of 
pixels from the MHD simulation is reduced by a factor of four 
such that the pixel scale is 23.8 km per pixel for both simulated 
snapshots. This is very close to the pixel size of the observations 
(see Sect. |2]) so we do not interpolate the synthetic data to a new 
scale. 

When comparing observations and simulations (see Sect.©, 
it is necessary to degrade the synthetic images to the quality of 
the observations. It is difl&cult to determine an exact degrada- 



tion PSF to apply to the synthetic images. The diflTraction by 
the finite telescope aperture as well as a low-order representa- 
tion of the phase aberrations are already included in the decon- 
volution performed in MOMFBD, and we deconvolve our ob- 
servational data using the diflferent theoretical interference filter 
PSFs as described in Sect. l3.2l Furthermore, some aspects of the 
image degradation (in particular stray light) are poorly known. 
However, the final reconstructed observations are noise filtered 
in a relatively aggressive way, eflTectively reducing the maximum 
achievable resolution from O'.'l to O'.'l? for the largest tilt angles 
(as described in Sect. l3.2l) . The fact that the noise filters are based 
solely on the power spectrum of the actual final reconstructed 
images means that they include unknown image-degradation ef- 
fects that are not compensated for in the observations. Therefore, 
for the intensity maps in Sect. O we use the tilt-angle-specific 
noise filters and, when calculating temperature gradients, we use 
two diflferent sets of noise filters (see Sect. 13. 2K as done in the 
observations. 

In addition to the above filtering, we add stray light to the 
synthe tic images in the sani e way and with the same parame- 
ters as lScharmer et al.l (|201 ih . using a r.'2-wide Gaussian stray- 
light PSF and considering a 58% stray light component at a wave- 
length of 538 nm. The large diflTerence in wavelengths between 
these two sets of observing data, in addition to the uncertain- 
ties already associated with this choice of straylight parameters, 
means that this is will not be the correct straylight function. 
However, our goal is only to make a morphological comparison 
between observations and synthetic data, and we do not attempt 
to remove straylight from the observations. The maps shown in 
Fig. [7] and Fig. [8] are computed from synthetic images degraded 
in this way (including straylight). 

The root mean square contrast (RMS) of the observed po- 
sition (7) quiet Sun is approximately that of the WB object 
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Table 3. RMS of the observational intensity versus the HD synthetics. 



Pos. 


Obs. 


Obs. 


HD 


HD 


(0 


MOMFBD 


filtered 




degraded 


(1) 


27% 


18% 


38% 


19% 


(2) 


18% 


13% 


34% 


16% 


(3) 


15% 


11% 


27% 


13% 


(4) 


12% 


10% 


20% 


11% 


(5) 


11% 


10% 


22% 


12% 


(6) 


11% 


10% 


25% 


13% 


(7) 


12% 


12% 


30% 


14% 



Notes. Columns in order: filter position, RMS of the quiet Sun patches 
shown in Fig.[T] RMS of the same patches after the most aggressive fil- 
tering (Sect.[3]2]), RMS of the synthesized HD images, and RMS of the 
degraded HD synthetics, including straylight as described in Sect. 14.41 

(which is 12.5%). After introduction of straylight, the RMS con- 
trast of the synthetic granulation images is still slightly higher 
than that of the observations (see Table [3]). The value of the de- 
graded near-continuum position (7) is consistent with compara- 
ble values fr om both synthetic and obs erved data discussed by 
I Wedemeyer -B ohm & Rouppe van der V oort ( 2009). 

We assume, for our aims, that the combination of this noise 
filtering and straylight model includes a reasonably realistic rep- 
resentation of the unknown image-degrading eflTects present in 
the observations. The fact that maps from both synthetic and ob- 
servational data (in Fig. [7] and Fig IS]) have the same dynamic 
range is a good indication of this and allows us to compare the 
morphology in the diflTerent computed quantities (see Sect©. 

5. Testing with a known 3D atmosphere 

To test our procedure, we use the undegraded synthetic im- 
ages produced from the field-free HD simulation as described 
in Sect. 14.41 We extract the temperature from the synthetic data 
using the methods in Sect. 14. II and 14. 21 and plot it against the ac- 
tual temperature from the input simulation snapshot, at the same 
column mass for each filter position in Fig. (6] Mean error and 
number of pixels below a certain error threshold are listed in 
Table IH To illustrate the procedure in detail, a plot of tempera- 
ture as a function of height with all the filter positions is shown 
for the coldest and hottest pixels in Fig.[5l 

There are no tilt angles or temperature ranges where the 
method fails entirely. The extracted temperature generally fol- 
lows the actual temperature from the simulation atmosphere. 

The errors are highest for the largest filter tilt angle (7) with 
an average over 300 K (see Tabled. The largest errors for this 
tilt angle correspond to pixels with very high vertical temper- 
ature gradient, where small errors in the depth determination 
will cause a big mismatch in temperature. Furthermore, this is 
the deepest position and thus has the highest overall temperature 
gradient. Uncertainties in the derived depth at this tilt angle in- 
clude the filter transmission profile being the broadest, Ca ii no 
longer being predominantly in the ground state (with the level 
population varying rapidly with temperature), and the contin- 
uum opacity being of significant importance. These last two as- 
pects are part of the limitation of the several one-point (height 
independent) assumptions in this simple analytical approach. 
Further one-point assumptions include the constant mean molec- 
ular weight and the one-point quadrature weight for the pressure 
integral (^i). We note that the continuum-opacity correction is 
computed for every filter position, but not as a function of tem- 
perature. The filter broadness introduces errors by mixing in- 



formation from diflTerent layers and by shifting the eflfectively 
sampled wavelength blue- ward of the fixed /lm,i in high-gradient 
areas and red- ward in low-gradient areas. The blue- ward shifting 
eflfect is expected to lead to an oflTset of temperature in the same 
direction of the offset observed in Fig.[6l 

The width of the filter also leads to a much narrower sampled 
height range than for spectroscopic data by causing the smallest 
filter tilt position (1) to have a greater contribution from the line 
wings than from the line core. The width also limits the number 
of relevant sampling points. In Fig. [71 chromospheric informa- 
tion is visibly mixed with the upper photosphere at position (1). 

The second worst position is position (5) with a mean error 
of 141 K. This position shows the same correlation of errors with 
high vertical temperature gradients as discussed previously, so 
the dominant error sources are probably the same as for position 
(7). The reason why position (6) is not worse than position (5) 
is that it is heavily contaminated with the All line (as visible 
in Fig.O, causing its formation height to be higher than what is 
calculated using Eq|4l This partially compensates the oflFset error 
from the depth-independent quantities. 

For all filter positions except (5) and (7), the mean error is 
below 100 K and more than half the pixels are below 100 K error. 
Somewhat surprisingly, for every filter position except position 
(5), the errors are smallest for the extremes of the temperature 
range. 

The MHD simulation is much colder than HolMul in the 
middle-to-upper photosphere. Due to the requirement of inten- 
sity calibration with the assumption that the observed quiet Sun 
regions will have the same average profile as the one computed 
from HolMul, testing this method with the MHD synthetic data 
leads to large oflfset errors and a generally large error scatter that 
may or may not be related to the intensity calibration. For this 
reason, we decided to focus on the HD test. Nonetheless, it is 
relevant to mention that generally, just like in the HD case, the 
error for the temperature extremes for each filter position is rela- 
tively small. This means that the temperature-extraction method 
seems to work well for pores and bright points. 

A similar test to the one shown here using simulation data 
and the Can K and H lines is performed by Sheminova (20^. 
She tested a manual fitting technique in full-resolution synthetic 
spectra and performed a test in selected columns from simulated 
3D atmospheres and another test in ID atmospheres averaged 
from those atmospheres. In the latter, she encountered problems 
mainly due to the non-linear temperature response of the Planck 
function, an issue that will always be present in unresolved struc- 
tures but is not very important at our high spatial resolution. In 
the resolved test, she encountered problems where the temper- 
ature gradient is sharp, similar to what we see. The fact that it 
is a manual technique precludes its use with filtergrams and, as 
the author points out, an automatic technique is desirable (also 
due to the sensitivity of a manual method to the initial assumed 
atmosphere). 

6. Comparison of observations and simulations 

6.1. Quiet Sun and HD synthetic data 

In Fig. [71 a relatively quiet area from the observations is com- 
pared with the HD simulation synthetic data for all the heights 
sampled. Since intensity and temperature maps look essentially 
the same, we show only intensity and temperature gradient maps. 
All the maps, at each height, have the same scale pairwise (each 
pair of synthetic data and observations). The intensity maps for 
the synthetics shown in the left column are not calibrated to 
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4000 5000 6000 7000 

Actual T [K]. Column mass density:2.17 [g cm"^] 



4000 5000 6000 7000 

Actual T [K]. Column mass denslty:2.60 [g cm"^] 



4000 5000 6000 7000 

Actual T [K]. Column mass density:3.03 [g cm"^] 




Fig. 6. Known simulated temperature from the 3D HD snapshot versus temperature extracted from the synthetic observations. The contour levels, 
from the center of the clouds to the periphery, correspond to 20%, 40%, 60%, and 80% of the number of points. 



Table 4. Results of the test in Sect.S 



Pos. 


< error > 


error < 50K 


error < \00K 


error < \00K 


(0 


(K) 






minus offset 


(1) 


73 


0.42 


0.73 


0.75 


(2) 


61 


0.51 


0.81 


0.81 


(3) 


60 


0.52 


0.82 


0.82 


(4) 


94 


0.30 


0.58 


0.76 


(5) 


141 


0.22 


0.38 


0.62 


(6) 


94 


0.29 


0.59 


0.71 


(7) 


323 


0.02 


0.04 


0.54 



Notes. The first column lists the position number as used in the text. The 
second column is the mean error from the test described in Sect. [5] and 
plotted in Fig.[6l The third and fourth columns are the fraction of pixels 
below 50K and lOOK of absolute error, respectively. The fifth column 
is the same as the fourth, but after a fixed temperature off'set has been 
subtracted. 



HolMul, so that the diflTerence in overall intensity is visible. In 
the right half of Fig. [71 temperature gradient maps computed as 



described in the previous sections are shown. The similarity be- 
tween observations and simulations is striking. 

The deepest layer (position (7) or 3.56 g cm"^) shows an ex- 
cellent match in both morphology and overall intensity level be- 
tween simulations and synthetics. However, what is remarkable 
is the similarity in the deep gradient maps (Gei at 3.30 g cm"^ 
and G57 at 3.09 g cm"^ ) as (7) is the filter position aflfected by the 
largest errors. Despite their large size, such errors are systematic 
and structure dependent, as seen in the previous section. Thus 
they aflfect structures present in both synthetics and observations 
in the same way. At this depth, both synthetic and observational 
temperature gradient maps show the interior of granules nearly 
uniformly as high-gradient structures and intergranular lanes as 
low-gradient structures. This is a pattern that holds up to the gra- 
dient maps at 1.15 g cm"^, with the diflTerence between granules 
and intergranular lanes gradually decreasing with height. 

The observational gradient map at position G56 is clearly 
noisier than the rest and appears very diflTerent from that of the 
simulations. Position (6) is polluted by the strong Al i blend, thus 
presenting mixed information from the different heights. The 
presence of the blend shifts the formation height higher, eflTec- 
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Synthetic Intensity Observational Intensity Synthetic Gradient Observational Gradient 



(1) 

0.17gcm-2 



(2) 

0.90 g cm-2 



(3) 

1.40gcm"2 





(4) 

2.18gcm-2 



(5) 
2.61 g 



(6) 

3.04gcm-2 



(7) 

3.56gcm-2 



I 1^1 ^ I I 1^ I I I I I I I I I 



0.54 g cm-2 



G23 
1.15gcm-2 



1.78 gcm- 




G45 
2.39 g cm-2 



G56 
2.82 g cm-2 



G51 
3.09 2cm-2 




Gei 
3.30gcm-2 



Fig. 7. Relatively quiet Sun patch (upper left corner of Fig. [T]) compared with the synthetic images from the HD snapshot. From left to right: 
intensity from synthetics, intensity from observations, gradient maps from synthetics and gradient maps from observations. The depths in terms of 
column mass shown on the left side are the average depth from all the pixels when computed as in Sect. 14.21 Each tick mark corresponds to one 
arcsec. Each pair of synthetic and observational images are scaled in the same way. 9 
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tively making it closer to position (5) than to position (7). This 
reduces the baseline from which gradients are computed, caus- 
ing noise to become more visible. 

Looking at the intensity maps for positions (6) to (4), at 
0.90 g cm"^ to 2.18 g cm"^, one sees that the transition from 
granulation to inverse granulation becomes apparent at a greater 
height in the simulations. 

Comparing the intensity maps in positions (2) and (1), we 
see that the pattern is nearly the same for both maps in the sim- 
ulation data. We note that the simulation does not aim at repro- 
ducing the chromosphere, so this is not surprising. However, for 
the observations, this pattern is also nearly the same but fuzzier 
for position (1). This is very suggestive of radiative scattering 
effects that are not captured by the 1.5 D radiative transfer used 
to produce the synthetic images. Scattering in the upp er pho- 
tospheric layers as observed in Can K was discussed by lRuttenl 
(1999) in the context of the observed fuzziness in network grains 
when compared with other diagnostics. While these layers are 
known to be highly dynamic, waves and heating mechanisms 
not present in the lower layers (such as the ones identified by 
ICarlsson eT al. (2010) and discussed in Rutten (2012)) seem to 
aff'ect mainly the shape of the intensity pattern. What we see in 
the area shown, however, is a constant pattern with an increase 
in "fuzziness". 

In the observational intensity maps, small bright points 
are visible that appear to have nearly the same size for all 
heights with the exception of the very diff'use line core position. 
Carefully inspecting the WB image, one can argue that the small 
bright points that are visible in the upper layers, including some 
as small as 5 pixels as measured from minimum to minimum, are 
about one pixel (0'.'034) wider than in the WB when looking at 
position (4) or higher in the atmosphere. This hints at expansion 
with height. However, this is close to our resolution limit, and it 
might be that these points are not resolved. 

Figure[7]illustrates that the highest possible spatial resolution 
is necessary to capture small-scale dynamics at these heights, 
even in the absence of magnetic fields. The morphological match 
between simulated data and observations in both intensity and 
temperature gradient maps strongly suggests that the physics of 
granulation is well captured by the HD simulation, but diff'er- 
ences exist between the synthetic images and observations when 
looking at the progression with height. 

6.2. An active area and MHD synthetic data 

In this section, we focus on interpreting Fig. [8l which shows a 
selected region of 10'' (~7.1 Mm) from the synthetic data pro- 
duced from the MHD simulation and a similar active area from 
the observations. This figure is built in the same way as Fig. [71 
However, due to space constraints and in order to show a larger 
image area that accommodates both pores and bright points, we 
cut out two rows. We choose to remove the images for position 
(6) and the gradients maps G^j and G56 due to the p ollution of 
position (6) from the Al i blend as stated in Sect. 16.11 The inten- 
sity map from position (1) is also removed because it is domi- 
nated by chromospheric features not present in the simulation. 
Both synthetics and observations present pores, and there are no 
point-like isolated bright features (like the ones visible in Fig. [7]) 
in either the observed or the synthetic data in the selected ar- 
eas. Even though the simulation atmosphere is much colder than 
the observations in the upper layers, which precludes detailed 
quantitative comparison, the morphological similarity between 
simulation and observations is impressive both in intensity and 
temperature gradient maps. 



For the deepest layer (G57 at 3.09 g cm"^) there is a region 
of strongly reduced temperature gradient at and surrounding the 
pores as well as bright points for both observations and synthet- 
ics (visible as a sprawling dark area in the gradient maps). This 
region is identifiable in the gradient maps up to 1 . 15 g cm"^ (G23 
), but between 2.39 g cm"^ and 1.15 g cm"^ (G45 and G23) small 
higher gradient structures (visible as bright filigree structures in 
the gradient maps) start rimming the outer border of this low- 
gradient area. These high-gradient structures expand and be- 
come stronger with height (relative to their surroundings). At the 
topmost gradient map (Gu at 0.54 g cm"^) for the synthetic data 
and at 1.15 g cm"^ (G23) for the observations, the expansion of 
these high-gradient structures leads to a near reversion of the pat- 
tern seen for the deepest layers (G57), with bright points showing 
higher temperature gradient than granules. Intergranular lanes 
are not distinguishable in the gradient maps at these heights (as 
was the case for the field-free HD simulation) for either synthet- 
ics or observations. The central part of the pores displays low 
gradient at all heights. 

In the deep layers, bright points have lower temperature gra- 
dient compared with their surroundings. A likely explanation for 
this is hot-wall radiation heating the mid-photo spheric optical 
depths inside fluxtubes, causing a "flatter gradient" as mentioned 
in [Ru tten ( 2010) when reproducing the hot- wall cartoon from 
[Spruiti(il976) . 

Nearl y all bright-point ch ains ("ribbons" and "flowers" as 
coined by [Berger et"aD ([20041) ) expand with height both in the 
synthetic and the observational intensity maps of Fig. [Sj This 
seems to be a difl'erent behavior than that seen in the isolated 
bright points mentioned in Sect. 16.11 The latter expand slightly, 
if at all, in the images, maybe because they are unresolved. 

A ringing artifact can be seen in the lower left part of the bot- 
tom gradient map for the synthetic data (at x = T\y = 3'' with 
origin in the lower left corner). This is purely due to the noise fil- 
tering and to the extreme sharpness of the feature. Diff'erentiating 
data enhances such artifacts. To a smaller degree, such ringing is 
also visible in the observed gradient maps but is hard to separate 
from the noise. The synthetic data are useful in drawing attention 
to such limitations of the data processing. 

The simulations replicate the observed chains of bright 
points and their expansion with height in intensity as well as 
pores. They also replicate the complex small-scale structure vis- 
ible in the observational temperature gradient maps. This degree 
of morphological similarity between simulated data and obser- 
vational data for both the MHD and HD cases makes it likely 
that both simulations capture the essential physics of the diff'er- 
ent observed structures. 

7. Conclusions 

Temperature stratification is extractable from filtergrams in the 
Can H wings at very high spatial resolution (up to O'.'IO) and 
over a wide atmospheric range in the same way as it has been 
used before in spectrograms. Temperature extraction was possi- 
ble with a simple analytical procedure making use of a filter with 
PSFs and passbands that vary with wavelength. 

The technique developed in this work to address non- 
simultaneity of frames in the context of image reconstruction 
using MOMFBD successfully minimizes spurious signal from 
the non- simultaneity of the recorded images and is already 
being adopted for other projects (e.g., 'Watanabe et al.[ ([2012[) : 
Sekse et al. (2012); Scharmer & Henriques ( 2012)). 

The limitations due to the broad and variable filter profile, as 
discussed in Sect. \5\ could be mitigated by introducing an addi- 
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tional rotation stage tiUing an interference filter further out into 
the Can H blue wing (for example, the 396.47 nm interference 
filter mentioned in Sect. O. This would minimize the tilt angle 
necessary to scan the whole wing. A more satisfactory solution 
would be to replace the setup with a tunable filter with narrow 
passbands and stable PSF (e.g. a Fabry-Perot). Significantly nar- 
rower passbands would also dramatically increase the sampled 
height range (see Sect. [5]) and allow proper separation of chro- 
mospheric and upper photospheric features that currently are 
mixed in the core filtergrams. 

A way to address the issues with the limitations of the 
temperature-extraction method used here could be to try to allow 
some quantities to be height dependent o r to iterate temperature 
in the deep layers as done for spectra by iRouppe van der Voori 
(12002"). However, in filtergrams, the range of the atmosphere to 
be iterated needs to be much larger. Thus, the ideal solution for 
this issue is to adapt a fully automati c inversion technique to th e 
Can H line wings such as t hose of | Bel lot Rub io et al.l (Il997h : 
ISocas-Navarro et al.1 ([T998I) ; Ide la Cruz Rodrigu ez et"al] (l2012h 
and lFrutiger et al.1 (|2000V 

The degree to which simulation data and observations match 
at several diff'erent heights in the Can H wing filtergrams is 
striking. This includes the morphology in diff'erential maps and 
the general progression of structures with height, excluding the 
smallest bright points. 

Temperature gradient maps are an interesting way to look at 
the morphology of structures and a testament to the success of 
the imaging-processing techniques presented in this work. The 
large diff'erence in gradient between areas with a strong mag- 
netic field and non-magnetic areas in the outer wings of Ca ii H 
makes these a possible proxy for automatically identifying both 
ma gnetic bright and dark features in a single map in deep lay- 
ers. 'Leenaartsetal] (l200d ) explored the blue wing of Can H as 
a proxy magnetometry diagnostic for bright points using diff'er- 
ential images from synthetic data at 396.73 nm and 396.48 nm 
(approximately the Ajn,i of position (2) and position (4) respec- 
tively). Our maps suggest that it is best to use diff'erentiation 
blue-ward of that range from the continuum u p to 3 96.43 nm 
when using this line as a proxy. Leenaarts et aD (!2006^ also con- 
cluded that Ha is the best proxy. However, the Ca ii H wing al- 
lows higher resolution. Additionally, it might be more interesting 
for instruments that target the entire Ca ii H line (for example, for 
its chromospheric potential plus context information) or that fo- 
cus on observing in the UV. Tha t would be the case, for example, 
of the upcoming IRIS mission (iLemen et al ] l20Tll) . which does 
not include polarimetry. 

In our upcoming work, we intend to interpret the temperature 
stratification obtained for the sunspot in the field of view. 

The coe fficient Fw = 1.33 x 10"^ is calculated by 
iRouppe van d er Voort ( 20021) from the broadening cross-section 
of lBarklem & O'Maral (il998h . 
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Fig. 8. An active patch with pores and ribbons of bright points (the pore region close to the vertex of the sunspot in Fig. [U compared with the 
synthetic images from the 3D MHD simulation snapshot from Lagerfjard & Nordlund (in prep.). From left to right: intensity from synthetics, 
intensity from observations, gradient maps from synthetics and gradient maps from observations. The depths in terms of column mass shown on 
the left side are the average depth from all the pixels when computed as in Sect. 14.21 Each tick-mark corresponds to one arcsec. Each pair of 
synthetic and observational images are scaled in the same way. 
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